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Abstract 

We survey the application of computer algebra in the context of 
gravitational theories. After some general remarks, we show of how 
to check the second Bianchi-identity by means of the Reduce pack- 
age Excalc. Subsequently we list some computer algebra systems and 
packages relevant to applications in gravitational physics. We con- 
clude by presenting a couple of typical examples. 
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1 Introduction 



Einstein's gravitational theory, general relativity (GR), is the valid theory 
for describing gravitational effects. In the search for making GR compati- 
ble with quantum theory and/or unifying it with the other interactions of 
nature (strong, electro- weak, superweak,. . . ), different schemes have been 
developed, like the gauge approach to gravity, including supergravity and 
metric-affme gravity, higher- dimensional Kaluza-Klein type models, string 
models, but also more conventional Hamiltonian (canonical) or Feynman 
quantization schemes or, more far-fetched, models based, e.g., on noncom- 
mutative spacetime geometries. 

The computer algebra (CA) programs applied in GR can be and partially 
have been extended to these more general frameworks. Still, it is probably 
true that most CA programs in gravity are applied in the context of GR 
followed by those for evaluating gravity-based Feynman integrals and for 
executing computations in the framework of gauge models encompassing non- 
Riemannian spacetimes. In this note we mainly concentrate on GR and will 
give a couple of examples. 

Detailed overviews of computer algebra in GR are given, e.g., by Brans @ , 
Hartley P|, Lake W7\, and MacCallum [ETJ, MM. 



2 Riemannian curvature in tensor and exte- 
rior calculus 

In GR and gravity, computer algebra was used as soon as it became available. 
The reason for this is that for solving standard problems it is required to 
manipulate a large number of terms and equations. We will clarify this by 
an example. A generic problem in gravity is to calculate the Ricci tensor from 
a given metric. A general form of a spacetime metric g in four dimensions is 
given by 10 independent functions = g^ of the coordinates ( 

3 

g = gij{x°,x 1 ,x 2 ,x 3 )dx i 0dx j . (1) 

i,j=0 
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The so-called Christoffel symbols are determined from the functions g^ by 
means of the following equations: 




. =r/. (2) 

m=0 \ / 

Here g km denotes the matrix reciprocal to gij. Since i,j, k are running from 
to 3, the Fij represent 64 functions. Because of the symmetry in i,j, only 



40 are independent. In the conventions of Schouten f31| , the Riemannian 



curvature tensor is derived from the Christoffel symbols in the following way: 

f) f ' f) r ' ^ 
p i _ _ 1 i fc _ _ 1 ik i \^ (r 1 r m — r 1 r m \ (i\ 

ni 3 k ~ Q x i Q x j V im i k L jm l-ik j ■ \o) 

Eventually we find the Ricci tensor as 

3 

Ric jk = Y,Rijk- (4) 
Now one can easily estimate that the number of terms in each of the com- 



ponents Ricjj may be very large. In [10] it is shown that in the general case 



this number is in the order of 10 000 for each of the components. Thus, only 
in simple cases can these calculations be done conveniently by hand. 

It would be most desirable to have a computer algebra system which 
allows to enter mathematical expressions in an analogous way as one would 
write them down on paper. That is, defining objects with abstract properties, 
doing calculations, and assigning explicit values to these objects should be 
possible in a natural way. Such kind of systems already exist. We illustrate 
this by an example. 

In terms of Cartan's calculus of exterior differential forms, Eq.(|3p can be 
displayed very compactly as 

Rj = dFj - r Q 7 A r/ . (5) 

Here, is the curvature 2-form and F^ the connection 1-form. The 
summation convention is assumed, i.e. summation is understood over the 
repeated index 7. Let us check the Bianchi identity 

dRj - IV A RyP + T/ A itV = (6) 
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on the computer. This can be deduced from Eq.(^]) and the properties of the 
exterior product and the exterior derivative. In Excalc, a Reduce package 
for exterior calculus, we first have to declare that all indices run from to 3. 
Then we declare T a ^ to be a 1-form and Rj 3 a 2-form. 

indexrange 0,1,2,3; 

pform gamma(a,b) = 1 , curv(a,b) = 2; 

Eq.(^|) and the left hand side of Eq.flf]) can almost literally be translated, 
curv(-a,b) := d gamma(-a,b) - gamma(-a,c) " gamma(-c,b); 

and 

d curv(-a,b) - gamma(-a,c) " curv(-c,b) + gamma(-c,b) ~ curv(-a,c) 

The negative sign in front of an index indicates that it is a subscript whereas 
a superscript is denoted by a positive (or no) sign. The last command yields 
zero, i.e. the Bianchi identity is confirmed. In [2§|, p. 234] it is shown how to 
handle this in the Mathematica package MathTensor. 



3 Abstract and component CA systems 

Packages (or systems) capable of performing symbolic calculations, as with 
the Bianchi identity above, are called abstract or indicial calculus systems Jl2| . 
They are necessary if one wants to investigate general properties of objects. 
So called component calculus systems are designed to calculate the compo- 
nents of unknown quantities from known ones. In our example, it was not 
necessary to introduce a basis or a metric in order to define and R a . 
In a typical component system, one first would have to enter the compo- 
nents of a metric. Then, by means of build-in routines, the components of 
the connection and curvature could be calculated. Both, abstract as well as 
component systems, allow to define new objects by means of various math- 
ematical operations like products or derivatives, e.g.. The difference here 
is that in the case of a component system it is always necessary to assume 
a specific basis and/or metric. In turn, there are abstract calculus systems 
which do not support computations of explicit values of components. Some 
packages, like MathTensor or Excalc, for instance, allow both, abstract and 
component calculations (see Table 1). In our Excalc example, we could as- 
sign explicit expressions to the components of T a ^. Then we would find the 
components of by calling curv(-a,b). 
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4 General versus special purpose systems 



Today it seems that most people use relativity packages of general purpose 
systems like Maple, Mathematica, Reduce, Macsyma and/or Derive. These 
programs offer very user-friendly front-ends and, moreover, a wealth of useful 
facilities, like simplification routines, programs for solving algebraic or dif- 
ferential equations exactly or numerically, Tex and Fortran interfaces, etc.. 

Special purpose systems are, as the name suggests, specialized to handle 
only a specific class of problems. Therefore the set of instructions is usu- 
ally very limited and the programming of these systems requires normally 
much more effort than is the case with general purpose systems. However, 
these systems are rather compact, very fast, and may sometimes be the only 
available facility to solve a problem. In the case of calculating Feynman dia- 
grams, e.g., special purpose systems like Schoonship or Form are often used, 
but there are also packages available for Mathematica and Reduce. A fairly 
widely used special purpose system for tensor calculus and general relativity 
is Sheep. For general relativity formulated in terms of orthonormal frames, 
the special purpose system Ortocartan has been recently updated [[[(J. 
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5 Applications 



Today many authors use computer algebra in order to obtain or confirm 
their results without mentioning it explicitly. In the following we present 
some articles which explicitly illustrate the applications of computer algebra 
in gravity. 

It was already stated that quantities like the Ricci tensor can reach an 
enormous size. Thus, for most applications (classification, numerics), these 
objects have to be put on the computer. The use of computer algebra has 
the advantage that there is no need to enter very large expressions for the 
curvature, e.g., but a comparably small input program which calculates these. 
Moreover, one can use the possibility to transmit programs or results by 
means of computer networks, cf . the closing remark in |7J . 

A standard application of computer algebra in GR is the classification of 
exact solutions. These are necessarily found in special coordinates. However, 
it is desirable to characterize the corresponding spacetimes in a coordinate 
independent way. This includes determination of Petrov and Segre types, 
calculations of curvature invariants (the trace of the Ricci tensor of Eq.(|]), 
for instance), the maximal isotropy group, etc.. The appropriate algorithms 
involve an enormous amount of calculational work. Programs for the widely 
used Petrov classification are available for most computer algebra systems. 
By means of the Sheep-package Classi |2l| , it was possible to create a search- 
able online-databank which includes nearly 200 exact solutions and their 
properties and which is still growing |[32| . 

To find out whether two solutions which look different are not just the 
same solution in different coordinates, one has to solve the so-called equiva- 
lence problem. This involves differentiation of the curvature tensor up to the 
seventh order. In |^TJ appropriate programs are presented for Sheep. 

Computer algebra is also very useful for finding new solutions of the Ein- 
stein equation. For vacuum, it reads Ricjj = 0. As we can recognize from 
Eq.(4) together with Eqs.(2) and (3), it represents a system of ten second 
order quasi-linear partial differential equations for the which are obvi- 
ously very difficult to solve. A simple example is given in || where the 
Schwarzschild solution is derived from a spherical symmetric line element 
by using the Reduce package Excalc. In it is illustrated how to use 
the Reduce-based system GRGec f° r searching for solutions of the Einstein- 



Maxwell equations. In [39] the Reduce package Classym was used to derive 



the Killing vector and Killing tensor equations and their integrability con- 
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ditions from a general form of a metric. Subsequently the Reduce package 
Crack has been used for solving these equations. 

In |34j] it is shown how to construct solutions of metric- affine gravity from 
solutions of the Einstein-Maxwell equations under heavy use of the computer 
algebra system Reduce. 

A further application is the derivation of a field equation from an action 
principle by means of variational calculations. A simple example is given 
in |28|, p. 303] where the Einstein vacuum equation is rederived. In [0, the 
gravitational field equation of a unified field theory of Pawlowski and Raczka 
were derived (and corrected) from the corresponding Lagrangian by means 
of MathTensor. 

In [f5?| an example is given of how to use the special purpose system Form 
to calculate Feynman diagrams in the context of quantum gravity. 

Another feature of many computer algebra systems is a Fortran interface 
which converts equations into Fortran readable form. This helps to develop 
programs for numeric calculations, as illustrated in |2|]. 

In || it is outlined of how to apply the computer algebra system GRTen- 
sorll to second-order black hole perturbations. 
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